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Understanding  Thermohaline 
Mixing  in  the  Agulhas  Return 
Current  from  Seismic  and 
Finestructure  Observations 

Ana  E.  Rice*,  Jeffrey  W.  Book  and  Warren  T.  Wood 


Naval  Research  Laboratory,  Stennis  Space  Center,  MS, 
United  States. 

*  Ana.Rice@nrlssc.navy.mil 


1  Introduction 

The  Agulhas  Current  System  in  the 
South  Indian  Oeean  (Figure  la)  is  one  of 
the  strongest  western  boundary  eurrent 
systems  in  the  world  due  to  its  large 
volume  transport  (order  70  Sv)  and  high 
distribution  of  mean  kinetie  energy  (up  to 
500  em^/s^  in  the  northern  Agulhas  Current) 
(Lutjeharms,  2007).  Downstream  of  the 
southwestward  flowing  Agulhas  Current, 
the  Agulhas  Return  Current  (ARC)  begins 
at  the  Agulhas  retrofleetion  in  the  southeast 
Atlantie  and  flows  eastward  in  a  quasi¬ 
stationary  meandering  pattern  between  38° 
and  40°  S  (Boebel  et  al,  2003).  Like  the 
Agulhas  Current,  the  ARC  exhibits  intense 
dynamies,  strong  mean  eurrents,  high  levels 
of  mesoseale  variability  and  abundant  eddy 
events  (Lutjeharms  and  Ansorge,  2001).  In 
addition,  the  elose  eontaet  between  sub¬ 
polar  and  sub-tropieal  water  masses  in  this 
region  produees  a  strong  frontal  zone  that  is 


eharaeterized  by  the  presenee  of  a  strong 
temperature  front  and  is  manifested  by 
interleaving  thermohaline  intrusions  at 
submesoseales. 

In  January-February  2012  the  Naval 
Researeh  Laboratory  (NRL)  and 
eollaborators  eondueted  a  field  experiment 
to  examine  frontal  zone  mixing  proeesses 
at  the  quasi-permanent  northward  meander 
of  the  ARC  in  the  vieinity  of  the  Agulhas 
Plateau  (Figure  lb).  Veloeity  results  from 
NRL’s  HYCOM  model  (Bleek,  2002) 
show  that  during  the  ARC  12  eruise 
(January  23-  February  8)  a  warm  eore 
antieyelone  that  had  been  loeated  over  the 
southeast  portion  of  the  Plateau  moved 
westward  and  briefly  merged  with  the 
ARC.  Both  seismie  and  traditional 
hydrographie  observations  were  eolleeted 
during  the  field  experiment.  Overall  8 
seismie  seetions  were  aequired  aeross  the 
ARC,  eddy,  and  zones  in  between  and 
these  were  eorroborated  with  simultaneous 
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XBT  casts.  The  seismic  observations  were 
complemented  by  repeating  the  transects 
with  18  CTD  stations,  37  casts  from  an 
Underway-CTD  system  (900  m  profiles 
taken  while  traveling  at  4  knots),  and  41 
micro  structure  casts  to  characterize 
turbulent  mixing  rates.  Vessel-mounted 
AD  CP  data  were  acquired  throughout  the 
seismic  and  hydrographic  surveys. 

This  study  focuses  on  analyzing  the 
seismic  and  hydrographic  observations  to 
elucidate  submesoscale  and  fine  scale 
frontal  zone  mixing  processes  by  assessing 
how  mesoscale  and  fine-scale  features 
influence  thermohaline  interleaving  and 
overall  mixing  processes  in  the  study 
region. 


(a) 


January  19,  2012 


February  8,  2012 


February  3,  2012 


Figure  1.  (a)  Sea  Surface  Height  (SSH)  in  the 
greater  Agulhas  System  on  January  19,  2012  as 
interpreted  from  satellite  data  using  NRL’s 
Modular  Ocean  Data  Analysis  System 
(MODAS)  (Fox  et  al.,  2002).  This  snapshot 
depicts  the  southward  flowing  Agulhas  Current, 
the  retroflection  and  Return  Current,  the 
presence  of  an  anticyclonic  eddy  on  the  Agulhas 


Plateau,  and  an  Agulhas  ring  west  of  the 
retroflection.  (b).  HYCOM  velocity  snapshots 
during  the  ARC  12  cruise  show  a  warm  core 
anticyclone  merging  with  the  ARC.  The  cruise 
track  and  location  of  seismic  sections  and 
hydrographic  stations  are  also  shown. 

2  The  seismic  oceanography  method 

The  seismic  oceanography  method  uses 
low-frequency  (10-200  Hz)  sound  waves  in 
the  ocean  water  column  to  generate  direct 
reflections  off  temperature  contrasts.  The 
generated  reflections  are  then  analyzed  to 
identify  water  mass  boundaries  of  fronts, 
eddies  and  oceanic  fine  structure  on  scales 
of  meters  to  tens  of  meters  (Ruddick  et  al., 
2009).  With  a  horizontal  and  vertical 
resolution  approaching  10  m  and  data 
spanning  from  about  50  m  depth  to  the 
seafloor,  the  technique  yields  band  limited 
measurements  of  ocean  reflectivity. 
Amplitude  peaks  and  troughs  in  the 
measured  seismogram  can  be  traced 
horizontally  for  many  10s  of  kilometers. 
Strong  temperature  gradients  produce 
strong  reflectivity,  which  when  convolved 
with  the  seismic  measurement  wavelet 
create  alternating  positive  and  negative 
values  in  the  seismic  image  from  single 
reflectors. 

3  Results  and  discussion 

During  the  ARC  12  field  experiment  the 
high  lateral  resolution  of  the  seismic 
observations  allowed  fine-scale  lateral 
tracking  of  thermal  intrusions  and  provided 
additional  information  about  the 
submesoscales  and  finestmcture  not 
discemable  by  the  hydrographic  data 
alone.  Figure  2  is  an  east-west  seismic 
image  collected  during  the  cruise  in  the 
section  that  crossed  both  the  northward 
flowing  ARC  and  the  southward  flowing 
portion  of  the  anticyclonic  eddy  (see 
Figure  lb).  The  vertical  temperature 
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Distance  (km) 


Location  Index  (CMP) 


Figure  2.  East-west  seismic  image  acquired  on  February  1,  2012.  In  this  section  the  northward 
flowing  ARC  in  the  vicinity  of  the  Agulhas  Plateau  interacted  with  the  southward  flowing  portion  of 
a  large  warm-core  anticyclonic  eddy. 


gradient  fields  derived  from  the  XBT’s  are 
also  plotted  for  eomparison.  Here  the  red 
(maximum  at  yellow)  and  blue  (maximum 
at  blaek)  indieate  positive  and  negative 
seismogram  amplitude  values, 

respeetively.  The  interaetion  between  the 
ARC  and  eddy  is  eaptured  in  these  aeoustie 
refleetions,  whieh  are  most  prominent  at 
the  edge  of  a  hyperbolie  flow  zone  (small 
region  of  weak  eurrents  in  between  two 
strong  eurrents  direeted  in  opposite 
direetions)  between  the  eddy  edge  and 
ARC.  There,  sub-polar  waters  interaet  with 
waters  of  sub-tropieal  origin  on  either  side. 
Refleetors  are  shallow  within  the 
hyperbolie  zone,  but  extend  for  many  tens 
of  km  away  from  the  zone,  reaehing  to 
depths  of  1100-1200  m  underneath  the 
ARC. 

Analysis  of  Temperature-Salinity  (TS) 
eurves  from  the  entire  hydrographie  dataset 
and  speeifieally  from  CTD  easts  taken  in 
the  ARC,  hyperbolie  zone  and  eddy  a  day 
prior  to  the  seismie  seetion  reveal  the 
nature  of  the  water  mass  boundaries 
assoeiated  with  the  observed  seismie 


refleetors  in  Figure  2.  The  water  eolumn  in 
the  ARC  features  warm  (<  25°C)  surfaee 
water  above  200  m  and  a  distinet 
subsurfaee  salinity  maximum  assoeiated 
with  South  Indian  Subtropieal  Water.  In 
the  frontal  region  between  the  ARC  and 
the  eddy  interleaving  is  evideneed  by  the 
presenee  of  strong  temperature  and  salinity 
gradients,  partieularly  between  200-800  m 
depth.  In  this  frontal  zone  eoneurrent  with 
the  hyperbolie  zone,  the  eharaeteristie 
salinity  maximum  at  200  m  no  longer 
exists  and  the  surfaee  water  within  the 
front  is  eolder  and  fresher  than  in  the  ARC 
(Figure  3).  A  spatial  map  of  salinities  from 
hydrographie  data  along  an  isopyenal 
reveals  that  this  fresh  water  may  be  from 
subpolar  origin  and  entrained  from  the 
southwest.  At  the  eddy  edge  TS 
eharaeteristies  are  a  eombination  of  TS 
eurves  in  the  ARC  and  TS  eurves  in  the 
frontal  region  and  indieate  that  the  eddy 
was  likely  spawned  from  the  ARC  at  an 
earlier  time  before  travelling  to  the  west 
and  interaeting  with  newer  ARC  waters. 
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Figure  3.  Temperature  and  salinity  profiles 
from  CTD  casts  in  the  ARC,  hyperbolic  flow 
zone  and  eddy  edge  taken  on  January  31,  2012 
along  the  seismic  transect  shown  in  Figure  2. 

An  analysis  of  a  separate  seismie 
survey  of  a  radial  erossing  of  the  eddy  (i.e. 
the  most  southern  east-west  seetion  on 
Figure  lb;  not  shown)  eombined  with 
simultaneous  XBT  profiles  eonfirms  a 
thermal  strueture  similar  to  the  ARC  and 
reveals  that  seismiely  measured  struetures 
underneath  the  eireulation  eore  and  outer 
edge  of  the  eddy  are  formed  by  thermal 
intrusions  that  ehange  in  shape,  size, 
strength  and  depth  over  10s  of  km  along  the 
seetion.  In  addition  seismiely  measured 
signatures  on  this  seetion  are  observed  to 
vary  aeeording  to  thiekness  of  the  intrusion 
due  to  the  eonvolution  of  the  seismie 
measurement  wavelet. 

Analysis  of  eonseeutive  XBT  easts  on 
the  seismie  eddy  radial  erossing  transeet 
(not  shown)  suggest  that  horizontal  traeking 
of  seismogram  amplitude  peaks  and  troughs 
do  generally  traee  eonneeted  thermal 
intrusion  features.  However  along  the 
seismie  transeet  that  erosses  the  ARC, 
hyperbolie  zone,  and  eddy  edge  (Figure  2), 
we  find  less  eonneetivity  between  thermal 
intrusion  struetures  in  eonseeutive  XBT 
easts,  partieularly  near  the  hyperbolie  zone. 
We  speeulate  that  vigorous  mixing 
dynamies  in  this  region,  partieularly  in  the 
loealized  region  of  the  hyperbolie  zone, 
eaused  thermal  intrusion  features  to  be 
mixed  away  over  horizontal  spatial  seales 


smaller  than  our  typieal  measurement 
spaeing  (5  km)  of  XBTs  there.  To  assess 
the  level  of  mixing  during  the  seismie 
survey  we  used  a  method  by  Holbrook  et  al. 
(2013)  whereby  turbulent  dissipation  of 
kinetie  energy  (e)  is  estimated  from  seismie 
horizontal  wavenumber  speetra.  However, 
the  signal-to-noise  ratio  of  our  data 
(aequired  by  a  medium-sized,  two  Gl-gun 
system),  was  insuffieient  for  this  analysis. 
Miero strueture  profiles  aequired  in  the 
hyperbolie  zone  at  a  different  time  however 
yield  estimates  of  s  as  high  as  10'^  W/kg'^ 
within  the  top  100  m  of  the  water  eolumn, 
whieh  suggests  the  presenee  of  aetive 
turbulenee  in  the  frontal  zone.  A  further 
horizontal  wavenumber  speetral  analysis  of 
the  vessel-mounted  ADCP  data  aequired 
during  the  seismie  seetion  reveals  the 
presenee  of  internal  waves  in  the  region  and 
suggests  that  the  elevated  turbulent 
dissipation  rates  may  be  due  to  these 
dynamies. 


Figure  4.  Temperature  and  salinity  from  CTD 
data  collected  during  the  ARC  12  cruise  plotted 
on  the  dT/dz  and  dT/dz  plane,  a  and  p  are 
thermal  and  haline  expansion  coefficients, 
respectively.  The  horizontal  dashed  lines 
indicate  rough  estimates  of  the  seismic  estimated 
detection  limits.  For  further  details  on  the  Turner 
angle  method  see  Ruddick  (1983). 

A  Turner  angle  analysis  employs  vertieal 
temperature  and  salinity  gradients  to 
quantify  the  relative  role  these  gradients 
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contribute  to  the  density  gradient.  Therefore 
Turner  angles  are  practical  indicators  of 
water  column  stability  to  double-diffusive 
activity  (Ruddick,  1983).  A  Turner  angle 
plot  (Figure  4)  computed  from 
hydrographic  data  can  also  be  used  for 
indicating  the  detectability  of  temperature 
gradients  by  the  seismic  method  in  a 
particular  study  area,  and  can  help  reveal 
the  nature  of  the  thermohaline  intrusions 
present.  The  seismic  method  is  sensitive  to 
vertical  temperature  gradients  in  the  water 
column,  but  the  wavelet  measurement 
convolution  combines  data  together  over  > 
40  m  in  the  vertical  and  therefore  small 
scale  gradients  in  temperature  that  are  weak 
become  indistinguishable  in  the  seismic 
measurements. 

For  our  seismic  observations  in  this 
study  we  estimate  a  rough  detection  limit  to 
be  of  the  order  of  0.04°C  per  meter.  The 
horizontal  dashed  lines  on  Figure  4  indicate 
these  rough  estimates  of  the  detection 
limits,  however  it  should  be  pointed  out 
that  all  data  above  or  below  both  sets  of 
lines  cannot  be  detected  by  the  seismic 
method  because  detection  is  a  function  of 
the  vertical  profile  of  dT/dz  in  addition  to 
the  strength  of  dT/dz  at  any  given  point.  In 
addition,  it  is  important  to  note  that  due  to 
the  thermal  gradient  detection  constraint, 
the  seismic  method  is  better  able  to  detect 
data  occurring  in  the  stable  (Turner  angles  - 
45°  to  45°),  and  salt  fingering  regimes 
(Turner  angles  45°  to  90°)  than  the  diffusive 
layering  regime  (Turner  angles  -45°  to  - 
90°).  In  the  ARC2  CTD  dataset  temperature 
gradients  >  0.04°C  per  meter  (see  data 
outside  the  dashed  lines  in  Figure  4)  data 
mostly  occur  in  the  salt  fingering  regime, 
with  a  smaller  subset  occurring  in  the 
diffusive  layering  regime,  indicating  that 
the  seismic  method  was  able  to  detect 


strong  temperature  gradients  from  these 
data.  We  speculate  the  gradients  conducive 
to  salt  fingering  in  the  study  area  were 
produced  as  a  result  of  warm  and  salty 
SISW  waters  existing  above  colder  and 
fresher  waters  below  in  the  ARC,  while 
diffusive  layering  gradients  were  produced 
as  a  result  of  interleaving. 
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